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ABSTRACT 

This work deals with the effect of reduced pressure 
in a controlled turbulent flow on electric discharge 
stabilization. A converging -diverging nozzle was designed 
and fabricated to provide a device for investigating this 
phenomenon. Turbulence was generated by perforated plates 
designed with a geometry that had proven to be successful 
in previous controlled turbulent work. Turbulence and 
discharge data are presented along with photographs of 
the discharge for no flow, laminar flow and turbulent flow. 
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I. INTRODUCTION 



Several classes of electric discharge gas laser devices 
are currently under investigation. (Ref. 1) A stagnant or 
no flow device accomplishes heat removal and discharge 
stabilization by diffusion to cavity walls. Maximum prac- 
tical power capability for this class of lasers is approxi- 
mately 500 watts (W) with a 10 meter (m) glass tube. Another 
one incorporates gas flow for convective cooling and dis- 
charge stabilization by diffusion to cavity walls. These 
devices have a practical power capability of 5kW with a 
7 m glass tube. A similar class, using convective cooling, 
increases discharge stabilization by various pre-ionization 
techniques with power capabilities greater than lOkW. 

Recent work with turbulent flow for discharge stabiliza- 
tion adds a fourth class of laser devices which have prom- 
ising power capability. B.iblarz and Nelson (Ref. 2) report 
significant discharge power increases with ambient pressure 
turbulent flow. Low pressure (80-200 Torr) turbulent flow 
studied by Wiegand and Nighan (Ref. 3) foimd increased 
discharge stabilization as a consequence of the turbulence. 
Another low pressure study by Wasserstrom, Crispin, 

Rom and Shwartz (Ref. 4) used a transverse secondary 
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flow to create vortical flow. They also found that vortex 
flow had a significant effect on discharge stabilization. 

The power capability of an electric discharge convec- 
tive laser is scalable with mass flow rate (Ref. 5) . The 
mass flow rate can be increased by increasing density, 
discharge area and velocity. Elevated pressure causes 
difficulty with discharge stabilization but more power can 
be coupled into the lasing mixture. No -flow breakdown 
data at different discharge pressures for 4 cm and 8 cm 
electrode gaps are presented by Loeb in Reference 6. Using 
these data, breakdown power for increasing pressure is 
shown in Figure 1. A single maximum is reached for a 
given gap length, indicating an optimum discharge pressure. 
The purpose of this project is to investigate the effect 
of sub -ambient controlled turbulence flow on discharge 
performance in the range above 450 Torr. 

Section II deals with selecting an economically fabri- 
cated device to explore the desired region and discusses 
the general philosophy supporting its design. A complete 
description of the entire apparatus is presented in Section 
III. The experiment was divided into two primary areas, 
turbulence tests and discharge tests, which are discussed 
in Section IV. The results of the experiment and conclusions. 
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Figure 1. Breakdown Power Versus Pressure (Loeb Data, 
Ref. 6) 
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with some recommendations for improvement, are found in 
Sections V and VI respectively. 

An interesting model of space charge convection effects 
is presented in Reference 7. Numerous experiments were 
conducted during this project to provide data to help 
establish and support the proposed model. Air was used 
for the discharge medium in all tests performed during this 
project. A CO^, and He mixture would be expected to 
produce superior results. 



12 



II. REDUCED PRESSURE SYSTEM DESIGN CRITERIA 



The design for reducing pressure in the electric dis- 
charge test section resulted from a preliminary study of 
two gasdynamic methods. First, the possibility of using 
a converging nozzle, such that the inlet was open to 
atmospheric air and the exit coupled to an evacuated 
reservoir, was considered. A pressure ratio of sixty per 
cent of atmospheric air could be obtained in the test 
section. However, the flow rate required to produce 
this pressure reduction combined with the need for at 
least a two minute run time dictated the use of a reservoir 
and vacuum pump much larger than was locally available. 

This approach, therefore, was subsequently abandoned. 

An alternative method which consists of using the flow 
characteristics of a converging-diverging nozzle was adopted. 
The system is driven by compressed air at the inlet, with 
the diffuser exit open to atmospheric air. A short, con- 
stant area section at the throat of the nozzle houses the 
electrodes for the electric discharge test section. The 
turbulence generating plates used in the experiment present 
approximately fifty per cent area blockage to the flow. 

To avoid choking the flow at the plate prior to reaching 
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desired test section conditions, a plate chamber was located 
upstream of the test section in the converging portion of 
the nozzle. 

The specific design of the nozzle was based on three 
required features. First, a pressure ratio of sixty per 
cent of atmospheric air and Mach 0.95 were the desired 
maximum operating conditions in the test section. Test 
section conditions would be reached when flow through the 
plate chamber, which housed a plate with fifty per cent 
area blockage, reaches Mach 0.3. The section linking the 
plate chamber with the test section had to be as short as 
possible, thereby keeping the turbulence generating plate 
within reasonable range of the test section. With this 
configuration discharge density would be a function of 
velocity but it was concluded that the effect of these 
two parameters on discharge stability could be separated. 
Details of the various sections are presented in Section III. 

It is known that vortex stretching occurs in accelerated 
flow of a wind-tunnel contraction, thus attenuating the low 
frequency, high energy turbulent spectrum of the flow. (Ref. 
8) The location of the turbulence generating plate inher- 
ently placed the produced turbulent flow in this adverse 
environment. Low frequency, high energy turbulence is 
known to have a stabilizing effect on the discharge. (Ref. 2) 
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It was felt that the plate hole dimensions could be altered 
to compensate for the attenuation and that the plate was 
located close enough to the test section to insure adequate 
turbulent flow through the discharge. This was later found 
to be incorrect. Turbulence decay was measured and is dis- 
cussed in Section V. 

Since breakdown conditions are slightly variable with 
changes in water vapor partial pressure (Ref. 9), a study 
was conducted to determine a method for drying the air 
supply. Chemical dryers provide the most effective means 
for drying moist air but commercial products capable of 
accommodating the expected flow rates are large, expensive 
units thus precluding their use. To provide consistency 
of measured data, the effect of relative humidity had to 
be taken into account. It was decided that breakdown 
voltage data would be corrected to a standard humidity 
of 0.6085 inches of mercury vapor pressure by applying a 
correction factor determined from curves presented in 
Reference 9. The procedure- used for determining the partial 
pressure of water vapor in the discharge region, and the 
correction factor curve, are presented in Appendix A. 

Recent work related to water vapor effects on breakdown 
voltage can be found in Reference 10. 
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III. EXPERIMENTAL APPARATUS 



The experimental apparatus is composed of three major 
subsystems. These subsystems include an air flow system, 
a variable high-voltage D.C. discharge circuit, and a 
turbulence measurement system. Figure 2 is a photograph 
of the laboratory and most of the equipment used in the 
experiment. 

A. THE FLOW SYSTEM 

The flow system consists of an air compressor, water- 
cooler heat exchanger, flow rate control valves, a plemmi 
section and a converging -diverging nozzle. Figure 3 is a 
schematic of the entire flow system. Air supply is pro- 
vided by a Carrier three-stage centrifugal compressor with 
a 4000 ft'^/min. capability at a maximiam pressure of two 
atmospheres. The flow is directed through a water-cooler 
heat exchanger that maintains a stabilized flow temperature 
at approximately 90^F. 

Three gate valves are used to control the flow rate 
through the test section. Two of these valves are located 
downstream of the heat exchanger and are the primary control 
valves for the flow. It was determined that the flow choked 
in the heat exchange when the primary valves were set for 
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Figure 2. Laboratory Equipment 
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Flow System Schematic 



max im um flow and this, therefore, necessitated the installa- 
tion of a heat exchanger by-pass controlled by the third 
gate valve. The by-pass eliminates the problem and full 
flow is available. 

A plexiglass plenum chamber used in Reference 11 was 
incorporated into this apparatus but it catastrophically 
failed during early operation and was subsequently replaced 
with a 3/4 in. thick, pl37Wood structure. The plywood 
plenum joints were glued and reinforced. Internal and 
external surfaces were sealed to prevent moisture absorp- 
tion. Aluminum honeycomb inserts dampen upstream turbulence 
and straighten the flow in the plenum. 

The converging-diverging nozzle is a plexiglass assembly 
of five individual sections. The assembled nozzle is shown 
in Figure 4. A converging section connects the original 
plexiglass plenum chamber to the plate chamber. The inside 
dimensions for the inlet, i.e., 9 in. x 9 in. were selected 
to integrate the new nozzle with the structure already 
present. The inside dimensions of the exit match that of 
the plate chamber. 

A plate chamber houses and supports the turbulence 

generating plates. The chamber's cross sectional area of 
2 

20 in. was selected so that at maximum operating conditions, 
i.e., Mach 0.3 flow in the plate section, the corresponding 
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Figure 4. Nozzle Assembly 






test section area (based on isentropic relations. Ref. 12) 
would be of sufficient size to house the electrodes. Re- 
moval of a pressure tight hatch allows easy access to the 
plate which, therefore, can be readily changed. The plate 
chamber is schematically illustrated in Figure 5. 

A short converging section links the plate chamber with 
the test section. Its length was determined by fabrication 
factors and the desire to have the plate as close as 
possible to the test section. The test section provides 
support for the electrodes and includes a number of features. 
The cathode can be moved to provide various electrode 
spacings, with a minimum spacing of 2.9 cm. Numerous 
threaded ports in the walls of the test section provide 
access for diagnostic probes which are secured in position 
by pressure tight fittings. The test section schematic is 
presented in Figure 6. A photograph of the fabricated 
section is shown in Figure 7. 

The diffuser section has an overall length of 41.18 
inches and a wall divergence angle of 2.5 degrees. Tufts 
were placed along the diffuser wall to check visually for 
flow separation. 
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Figure 7. Test Section 
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B. THE DISCHARGE CIRCUIT 



Three primary elements comprised the variable voltage, 
D.C. discharge circuit. These elements include a variable 
high-voltage, D.C. source, the discharge electrodes and 
current/voltage meters. Figure 8 is a schematic of the 
circuit. 

Two high-voltage power supplies were available and 
were independently used in the experiment. A Sorenson 
high-voltage, D.C. power supply with a maximum output of 
30 kilovolts (kV) and 20 milliamperes (mA) powered the 
electric discharge in most of the tests. An old Sorenson 
Beta high-voltage, D.C. power supply with a rated maximum 
output of 60kV and 50mA was also used but found capable 
of producing only 34kV and 18mA. The power supplies were 
interchanged to cross check measurements of experimental 
discharge power data. Excellent agreement was observed. 

The discharge electrodes consist of a pin anode arrange- 
ment similar to Aunchman's (Ref. 11) and a wire grid cathode. 
The anode pins were made from 0.0625 inch steel wire and 
mounted on an airfoil machined from 0.5 to 0.125 inch brass 
strip. Figure 9 shows the details of the pin electrode. 

The wire grid cathode arrangement resulted from a series 
of preliminary tests utilizing Aunchman*s equipment. The 
airfoil cathode used by Aunchman was replaced by an aluminum 
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Figure 8. High Voltage Discharge Circuit 





Figure 9. Pin Anode Schematic 
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honeycomb cathode to establish a more uniform field, and 
did perform slightly better than the airfoil arrangement. 

The honeycomb cathode was not rigid enough to withstand 
the flow forces and hence was replaced by a wire grid 
cathode which performed just as well. It was learned later 
that wire grid cathodes are the preferred cathode geometries 
of axial flow systems. (Ref. 13) A photograph of the elec- 
trodes is shown in Figure 10. Mounted on the pin anode is 
plate III used in the tests and described later. 

Current measurements were taken with three Triplett 
milliammeters with ranges from 0.0 to 1.0, 10.0 and 50.0 
milliamperes wired in series through a "make before break" 
switch. A Simpson 260 V-O-M was also used for a cross 
check measurement of current. Originally, discharge voltage 
was measured with a Sensitive Research electrostatic volt- 
meter but a short circuit was detected in the meter and 
precluded its use. The power supply control panel voltmeter 
readings are within one per cent of the Sensitive Research 
meter and it was therefore used to measure the discharge 
voltage. This error was taken into account when an error 
bar for measured breakdown voltage was determined. High- 
voltage leads were used throughout the circuit. Power 
supply connections were made through a highly polished, oil 
submerged, brass sphere. 
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Figure 10. Electrodes 
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C. TURBULENCE EQUIPMENT 

Turbulence diagnostic equipment includes a Thermo 
Systems Inc. hot wire anemometer, a Hewlett-Packard oscillo- 
scope, a General Radio Company, Type 1921, real time 
analyzer and a Hewlett-Packard 7035B X-Y recorder. Hot 
wire probe filaments are 0.00015 inch diameter tungsten 
wires spot welded to the probe supporting pins. Hot wire 
anemometer components include a model 1051-2 monitor and 
power supply and a model 1054A-20 anemometer module. 

The oscilloscope was used for setting and adjusting 
the stability and frequency response of the anemometer. 
Difficulty was encountered in matching probe /anemometer 
impedance following the procedure outlined in the operators 
manual. To alleviate the problem, various lengths of co- 
axial cables were kept on hand for coarse impedance matching, 
followed by fine adjustments of the anemometer. In general, 
the hot-wire system proved difficult to use. Considerable 
skill and time were required to spot weld the tungsten 
wires in place. The wired probes were fragile and required 
frequent replacement. 

The real time analyzer is composed of a model 1925, 

30 channel multifilter and a model 1926 multichannel RMS 
detector. Based on experience gained in Reference 2, 
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four-second integration times were used during all measure- 
ments, The X-Y recorder enables rapid and accurate data 
collection from the memory output of the analyzer. All 
equipment is grounded through a separate grounding system 
in the laboratory. 

D. OTHER DIAGNOSTIC EQUIPMENT 

A number of diagnostic instruments used in the experi- 
ment were fabricated locally. Iron-constantan and copper- 
constantan thermocouples were used for total temperature 
measurements. A pitot-static probe, using either water 
or mercury manometers, provided pressure data. The thermo- 
couples and pitot-static probe were carefully calibrated 
at a local facility. Pressure and temperature measurements 
of a free jet from Mach 0 to Mach 1.0 were compared with 
highly accurate reference probes. The calibration pro- 
cedure resulted in detection of a faulty probe, which was 
discarded, and permitted a high level of confidence in 
experimental measurements. Other diagnostic equipment was 
also calibrated according to operating manuals. Calibra- 
tion results were used to estimate an error bar for the 
various measured test parameters. Error analysis is 
presented in Appendix B. 
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IV. EXPERIMENTS CONDUCTED 



The experiments performed were divided into two main 
areas. First, turbulence spectrum sampling was required to 
investigate the amount of turbulence decay encountered with 
accelerated flow in this apparatus and to make necessary 
modifications. Large eddies must be present in the dis- 
charge region to have a stabilizing effect in breakdown 
as evidenced by Reference 1 and Reference 4. Second, maxi- 
mum discharge performance was to be investigated with 
turbulent and non -turbulent flow. The diffuser was removed 
to obtain reference data at near atmospheric conditions. 
With the diffuser on, density reduction effects could then 
be studied. 

A. TURBULENCE TESTS 

In Reference 11, spectrum distributions were obtained 
for various plates. Plate 9 had produced the most intense, 
homogeneous turbulence distribution. However, discharge 
performance tests at a velocity of 200 ft/sec. (61 m/sec.) 
showed Plate 1-6 outperformed Plate 9. Preliminary tests 
using these plates and Aunchman's test section, showed 
that Plate 9 was superior at higher velocities. Plate 9 
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geometry was adopted for use in the current apparatus. 

The various plate modifications used are identified as 
Plates I, II and III. Figure 11 describes these plates 
in detail. 

A spectrum distribution of the turbulence generated 
by the pin anode was recorded at a test section velocity 
of 200 ft/sec. (61 m/sec.). Plate I was used to investi- 
gate spectral decay associated with flow acceleration. 

This was accomplished by placing Plate I at the upstream 
end of the plate chamber. (Position A, Fig. 3) . Turbulence 
spectrum was measured in the chamber 3.5 inches (8.89 cm) 
downstream of the plate. Plate I was then relocated at 
the downstream end of the chamber (Position B, Fig. 3) and 
turbulence in the test section 3.5 inches (8.89 cm) down- 
stream of the plate was recorded. Comparing the spectra 
obtained by this procedure, acceleration effects on turbu- 
lence were determined. Spectrum decay was much greater 
than initially anticipated. Plate II was manufactured and 
tested in the same manner. With Plate II installed at 
location B, spectra in the test section were measured with 
and without the diffuser. 
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Figure 11. Turbulence Generating Plates 
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B. DISCHARGE PERTORMANCE TESTS 



This part of the experiment involved discharge power 
performance under the influence of the various plate and 
nozzle arrangements used in the turbulence tests. Three 
internal configurations were studied with the diffuser 
removed. First, maximum discharge performance with in- 
creasing velocity was measured for a no-plate, "laminar" 
condition. Again, maximum discharge performance with in- 
creasing velocity was recorded with Plate 11 installed in 
the plate chamber. Since power supply limitations were 
reached prior to obtaining optimum desired test section 
conditions, it was decided to manufacture Plate III and 
mount it on the immediate upstream edge of the pin anode. 
This configuration is shown in Figure 10 and the plate 
position is identified with position C in Figure 3. It 
was anticipated that the flow would choke at the plate 
prior to optimum discharge conditions, but it was felt that 
power supply limitations would be reached prior to choking. 
Attempts were made to obtain a turbulence spectrum of Plate 
III in the discharge region, but hot wire tungsten fila- 
ments could not survive this intense turbulent flow even 
at low velocities. All the above tests were repeated with 
the diffuser on. 
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A number of other tests involving variable voltage 
rates and various gap lengths were performed to provide 
a data basis for a model presented in Reference 7. 
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V. DISCUSSION OF RESULTS 



Spectral data indicate that a severe low frequency 
turbulence decay occurs for plates positioned at location 
B (Fig. 3) in the accelerated flow. These data are pre- 
sented in Figure 12. Comparing the unaccelerated spectrum 
with the accelerated spectrum of Plate I illustrates the 
decay encountered. Spectra shifted slightly to the right 
with increasing velocity as was observed in Reference 11. 
Plate III turbulence spectrum could not be measured, for 
reasons previously stated, but some partial data did show 
turbulence intensities above 40dB. No changes in the 
spectrum were observed with diffuser on and diffuser off 
nozzle configurations at a test section velocity of 
200 ft/sec. (61 m/sec.). 

Total mass flow rate through the nozzle was computed 
from plate chamber diagnostics. Mass flow through the 
discharge region was less than total mass flow because of 
plate and pin anode blockage. The plate does not extend 
from wall to wall. The width of the plate is less than 
the total width of the test section and thus provides a 
"boundary layer" slot. The slot allows higher mass flow 
along the walls to remove discharge waste heat and to 
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PLATE CHAMBER 
O Plate I, V = 100 ft/sec 

TEST SECTION 

■a Plate I, V = 179 ft/sec 
<$► Plate II, V = 200 ft/sec 




FREQUENCY (Hz) 

Figure 12. Turbulence Spectra 
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remove any water condensation. Mass flow outside the dis- 
charge volume was determined from diagnostics in that 
region of the test section and subtracted from the total 
mass flow giving the mass flow rate through the discharge 
volume. Figure 13 was constructed from mass flow rate 
data. 

An error analysis for discharge performance curves is 
presented in Appendix B. The magnitude of the computed 
error bar was extremely sensitive to pressure readings. 

At high pressures, mercury manometers were used producing 
a larger error bar than the readings from the water mano- 
meters . 

Maximum discharge power per mass flow rate or energy 
per mass obtained for no plate "laminar" flow and turbulent 
flow using Plate III located at position C (Fig. 3) are 
presented in Figure 14. These curves show an initial de- 
crease in energy per mass to a minimim and then rapidly 
increase. At low mass flow rates the total power remained 
relatively constant causing the initial decrease. As mass 
flow rate increased, convective cooling, turbulence and 
sub-ambient density increased the stability of the dis- 
charge resulting in greater total discharge power. The 
model presented in Reference 7 suggests how velocity. 
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Figure 13. Mass Flow Rate Through the Discharge Versus Total Mass Flow Rate 



POWER PER MASS FLOW RATE (kW/kg/sec) 




Figure 14. Power Per Mass Flow Rate Versus Mass Flow 
Rate Through the Discharge 
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turbulence and density influence discharge stability. 

The data presented here help support the validity of the 
proposed model. 

Sub -ambient density effects on discharge stability 
caused an increase in energy per mass as was anticipated. 
The advantage of turbulent, sub-ambient flow can be 
realized from these curves. At the power supply current 
limited points, i.e., the last data point for each curve, 
more than twice the energy per unit mass was obtained for 
sub -ambient turbulent flow at about half the mass flow 
rate when compared with atmospheric "laminar" flow. 
Another way of looking at this is to superimpose lines of 
constant power on the curves of Figure 14 as shown in 
Figure 15. A given power level can be reached with sub- 
ambient turbulent flow at a much lower "cost" in mass 
flow rate. 

Figure 14 also shows that the laminar curves appear 
to be approaching an upper bound on energy per unit mass 
while the turbulent curves are rapidly rising. Data from 
plates I and II were identical to the laminar curves, 
demonstrating the severe turbulence attenuation of the 
accelerated flow. 
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POWER PER MASS FLOW RATE (kW/kg/sec) 




Figure 15. Constant Power Curves 
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An average discharge velocity was computed from dis- 
charge mass flow rates. Convective cooling and turbulent 
effects on discharge stability with increasing average 
discharge velocity can be seen in Figure 16. Increasing 
velocity provided convective cooling which increased dis- 
charge stability and hence greater power density. Higher 
power densities were obtained with turbulent flow at 
velocities greater than 60 m/sec. The slope of the 
"laminar" curve appears to be decreasing above 160 m/sec. 
Peak power would probably be reached at choked conditions. 
The extreme ends of the curves are power supply current 
limits, not energy density limits. Figures 17 and 18 
give breakdown voltage and breakdown current data for 
increasing discharge velocity. 

Photographs of the discharge near breakdown conditions 
taken for no -flow, "laminar" flow and turbulent flow are 
presented in Figure 19. The no-flow photograph was a ten- 
second exposure. The general appearance of the glow 
discharge shows a low intensity discharge extending be- 
tween each pin and the cathode. The bright region at the 
bottom of the photograph resulted from a breakdown spark 
during the time exposure. In contrast, the "laminar" and 
turbulent flow photographs were one-second exposures. In 
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POWER DENSITY (W/cm-^) 



• Diffuser off. No Plate 

★ Diffuser on. No Plate 



O Diffuser off, Plate III 




Figure 16. Power Density Versus Average Velocity 
Through the Discharge 
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Figure 17. Breakdown Voltage Versus Average 
Velocity Through the Discharge 
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ANGLE OF VIEW NO -FLOW LAMINAR FLOW TURBULENT FLOW 





the "laminar" flow photograph, the individual pin dis- 
charges were more diffused. An intense, short, narrow 
glow extending from the pin anode tips had developed. 

This is a region of high current density. As shown in 
the turbulent photograph, the discharge had become even 
more diffused. The bright region at the tips of each pin 
was wider and visibly more dispersed. The localized 
current density in this region has become more diffused 
due to turbulent mixing. Similar observations are 
reported in Reference 3. The outer or periphery pins of 
the anode were the first pins to develop the intense bright 
region described above. The center pins did not "light 
up" as brightly as the outer pins. 
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VI. CONCLUSIONS AND RECOMMENDATIONS 



Even though maximum operating conditions were not 
reached due to power supply current limitations, sub- 
ambient, controlled turbulent flow had a significant 
effect on discharge stability in the working range of 
this experiment. Turbulence generating plates must be 
kept in the immediate vicinity of the discharge (location 
C, Fig. 3) to avoid spectral decay in a wind-tunnel con- 
traction. Power input is scalable with mass flow rate. 

The combined effect of reduced pressure and turbulent flow 
at maximum performance obtained in this experiment was to 
increase discharge power by 750% over ambient laminar 
flow at the same mass flow rate. 

The upper limit of power input for this apparatus 
would occur under choke conditions with Plate III installed 
at location C, but a power supply with sufficient output 
would be needed to reach that level. The slope of the 
turbulent flow curves may or may not eventually decrease 
in the same manner as the "laminar" curves, and this 
needs to be further investigated. Turbulence spectrum 
data could be enhanced by laser velocimeter measurements 
of the turbulent flow. 
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At a given average discharge velocity the total dis- 
charge power can be obtained from Figure 16 for either 
the turbulent or "laminar" case. This power is represented 
by a line of constant power in Figure 15. The intercepts 
of this constant power line give mass flow rates correspond- 
ing to the diffuser on or off configuration. Since dis- 
charge area and velocities are equal, only density accounts 
for the difference in mass flow rates. This conclusion 
is supported by the data. 

By extrapolating the probable path of the "laminar" 
curves in Figure 14 (see Fig. 20), decreasing density ratio 
versus total discharge power can be projected. The ratio 
of mass flow rates corresponding to the intercepts of 
constant power lines gives p / . The resulting curve is 

shown in Figure 21. 

There are a number of improvements that can be incor- 
porated into this apparatus. As mentioned earlier, a 
larger power supply would vastly improve the capability of 
the system. This would also require additional electric 
shielding of the high voltage circuit to protect the 
operators. The center row of anode pins had a glow inten- 
sity less than the periphery pins for all flow and no -flow 
tests. This may have been influenced by the cathode wire 
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Figure 20. Extrapolated Laminar Curves 



■m 



POWER (kW) 



• Extrapolated Points from Fig. 20 
O Experimental Data 

1 , 0 -- 

. 9 -- 






. 7 - 

. 6 -* 



. 5 -- 



. 4 -- 

. 3 -' 



. 2 -- 



. 1 -- 



0 ^ 

.6 .7 






# 



# 






o 



o 



o 

o 



.8 



- 4 - 

.9 



-9 

1.0 



p. 



Figure 21. Power Versus Density Ratio 
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grid frame. Optimizing pin density for turbulent flow may 
also produce interesting results. An electrode gap length 
of 2,9 cm was used for all discharge data presented. With 
a larger power supply the optimum gap for the corresponding 
operating density of the discharge region could be experi- 
mentally determined and improve the performance. 

A closed, pressure tight system as suggested by Reference 
11 would alleviate the problem of relative hxamidity and would 
provide independent control of density and velocity. A 
closed system would also enable economic use of valuable 
lasing gas mixtures. 
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APPENDIX A 



PARTIAL PRESSURE OF WATER VAPOR IN DISCHARGE 



The relative hiomidity and temperature of ambient air 
determines the partial pressure of water vapor at the 
compressor inlet by using 
( 1 ) 

where P^ is the partial pressure of water vapor, (D is the 
relative humidity and Pg is the saturation pressure deter- 
mined from steam tables. A detailed steam table can be 
found in Reference 14. From Dalton's Law, 

(2) P = + Pa 

where the subscript a refers to the partial pressure of 
dry air. The specific humidity is given by 



(3) 



r 



£w ^ _ 

Pa 



0.621 



LB OF WATER 
LB OF AIR 



and are the molecular weights of water and air 
respectively. 

From Eqs. (1) and (2), and measured ambient pressure, 

P^ can be determined, which in turn, when substituted into 
Eq. (3) will give the specific humidity of the flow. The 
specific humidity will remain unchanged provided condensation 
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does not occur. If the flow temperature drops below the 
local dew point, condensation will take place and, therefore 
reduce V . 

To determine the partial pressure of water vapor in the 
discharge region, the static pressure and temperature of 
the test section were measured. Assuming for the present 
r remained constant, Eqs. (2) and (3) are used to deter- 
mine in the test section. 



(4) 



P 

(1 + r /0.621) 



(5) 



Pr. = 



0.621 



The steam tables were used to determine P^ in the test 

O 

section and substituted into Eq. (1) to determine the rela- 
tive humidity. If the calculated ® was greater than or 
equal to one, F had decreased or was about to and the 
partial pressure of water vapor in the discharge region, 

P^, was equal to Pg. If ® was less than 1, the partial 
pressure of water vapor as calculated in Eq. (5) was correct 
A corresponding breakdown correction factor was then 
determined from Figure 22. All data was corrected to 
standard humidity conditions. 
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APPENDIX B 



ERROR ANALYSIS (REF. 15) 

The uncertainty of the results depended on accurate 
measurement readings of breakdown voltage and current, 
total flow temperature and static and stagnation pressures. 
The following is an estimation of the uncertainty interval 
of these variables. 

Breakdown voltage: 

Vb = (V ± 0.2) 0.99 kV 

The 0.99 comes from using the control panel 
voltmeter. Breakdown voltage was first 
corrected to standard relative humidity. 

Breakdown current: 

lb == 1 ± 0.2 mA for 0 < I < 10 mA 

I, = I + 1.0 mA for 10 < I < 24 mA 

D — 

Total Temperature: 

Tq = T + 2 °R 

Stagnation Pressure: 

= p + 0.005 PSIA for H 2 O manometers 
P = p + 0.02 PSIA for Hg manometers 
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Static Pressure: 



Pg = p + 0.005 PSIA for H 2 O manometers 
Pg = p + 0.02 PSIA for Hg manometers 

The intermediate results are fimctions of the test variables 
as stated below. 

Power : 

P = VI watts 



Flow Temperature: 

T = T(Tq, Pq, Pg) isentropic 

Velocity: 

V = v(Tp, Pq, Pg) isentropic 



Density: 

P =P(Tq, Pq, Pg) isentropic 

First, consider the vincertainty in experimental values of 
power. 



1) P = VI watts 

2) P + aP = 0.99 (V + A V) (I + A I) 
Subtracting eq. 1) from 2) gives 

^P = 0.99 (VI + V AI + I aV+ aV aI)-1 
A P = -0.01 VI + 0.99 (V A I + I A V) 



A P 



= 0.99 



,Ai + AV) 
^ I ^ V 



0.01 
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In a similar manner the uncertainty interval was 
determined for the remaining test quantities and are 



summarized below. 



Flow Temperature: 

To + ATo 

AT 

— = / Po+ I aP„I 

\ Ps± I iPs' 



Velocity: 

^ J j 12015.08 (T„+ A To) 1^1 

V ^ft/secj 
Density: 



A P 


2.70 


’ Ps ± 1 


APsI 


1 -1 


ribml 

P lf^\ 


P 


~ P 


_ T + 1 


A T 1 


J ’ 



Mass Flow Rate: 

AM A p Ay w r Ibm l 

M * P V ’ LsecJ 



0.5 

/ > A Ps* \ 1 I -1 

VPo± I A PqI / J ) 
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